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SUBROUTINE  MXPD1NPP1  , N«fH  , ALPHA , BETA , SIOSQ, NOBS  , Y , BEL  , NTF  , 

1 NTL , NVF , N VL , NN , I ROUS 1 , 1 R0US2 , AL  , D , UKL , R , RO , E , X , YPD , I FAULT  > 

THIS  SUBROUTINE  CALCULATES  LEAST  SQUARES  PREDICTORS  FOR  A NIXED 
AUTOREGRESSIVE  NOVINO  AVERAGE  PROCESS  OF  ORDER  <NP,NB> 


I 


DIMENSION  ALPMA(NPP1 >,BETA(N0P1 >, Y < NOBS >, AL < IROUSI , IROUS1 >, 
1 D< IROUS1 ) . UKL < NQP1 ) , R (NflPI > , E ( IROUS2  > . X ( NOBS  > , YPD(NN) 

DATA  ZERO, IOPTB, IOPTE/O. 0 , 0, 1/ 

t 

TEST  FOR  INVALID  PARAMETERS 
IFAULT-1 

IF(NPPt .LT.1 > 00  TO  16B 
IFAULT-2 

IFtNQPt .LT.1 > 00  TO  160 
IFAULT-.1 

IF<<NTL-NTF*1)»(NVL-NVF*1).0T.NN)  00  TO  160 
I FAULT -6 

IFC IROUSI. LT.N0P1)  00  TO  160 

FIND  E ARRAY 

NO-NQPI-i 

NP-NPP1-1 

NPPNO*  NP  +NQ 

ir(NO.FO.O)  00  TO  40 

NX-NOBS-NP 

DO  20  1-1. NX 

NP1'I-NP+I 

C-YCNFPI > 

DO  10  .1-1  ,NP 

NI’PIN.I'NPPI-J 

r.  r<AI  FMAC  J>*Y<Nf'PIH.I> 

10  MINI  I HUE 

xu>-r 

20  rONTINI* 


. 


CA1  L HAORTMI  IUPTB.  IOPTE  , NO,  DEL  , NX  , NX  , X , f<ff  TA,  SIGSQ,  I RUNS 1 , 
1R.RO.NM  , D.AL  ,E,IF1  ) 

C 

IFAULT-5 

IF < IF1 . EQ.  5 > GO  TO  160 
DO  30  1*1 , NX 
N1-N0BS-I+1 
N2*N1-NP 
F < N1 >-E(N2> 

30  CONTINUE 
C 

C FIND  PREDICTORS 
C 

40  CONTINUE 

N1-NVL-NVF+1 

N2*NP*NVF 

DO  150  NT-NTF , NTL 

N3*<NT-NTF>»N1 

NTEHP 1 *NT -NP-NQ 

NTENP*NT-NP 

IF ( < NF  PNQ.EQ. 0/ -OR. (NT.LT. 1 > .OR. (NT. GT.NOBS) ) GO  TO  55 
IF(,<NP.NE.OKAND.  < NT . LT . NPPNQ) > GO  TO  55 
GO  TO  60 
55  CONTINUE 

DO  50  1-1 ,N1 
N3PI*N3*I 
YPD(N3PI)-ZER0 
50  CONTINUE 

GO  TO  150 
60  CONTINUE 

IF(NP.EQ.O)  GO  TO  80 
DO  70  1-1 , NP 
N4-NTEHP+I 
X< I >-Y  <N4)  - 
70  CONTINUE 


80  CONTINUE 

DO  130  NV-1,NVL  • 

C-ZERO 

IF( (NQ.EQ.O) .OR. <NV.OT.NQ>>  GO  TO  100 
NTPNV-NT+NV 

IF< (NP.EQ.O) .AND. < NTPNV. LE . NQ) ) GO  TO  100 
DO  70  K-NV, NQ 
KK-NTPNV-K 
C-C ♦BE  TA(K)»E(KK) 

70  CONTINUE 

100  CONTINUE 

NPPNV-NP+NV 
IF(hP.EO.O)  GO  TO  120 
DO  110  1-1. NP 
NP1-NPPNV-I 
C*C-Al_PHA<  I )*X  < NP1  ) 

110  CONTINUE 

120  CONTINUE 

X ( NPPNV) -C 
130  CONTINUE 

DO  140  NV* 1 , N1 

N3PNV*N3*NV 

N5*N?«NV-1 

YPD  < N3PNV ) * X < N5 ) 

140  CONTINUE 

1^  CONTINUE 
C 

irAUtT*0 
160  CONTINUE 
RETURN 
END 

SUBROUTINE  NAORTM( IOPTB , TOPTF , NQ , DEI  , ITERS , IR0US1 . Y , BETA, 

1SIU.VU,  IRONS?  , R , RO  , MKI. , D , AL  , E , IF  AUL  T ) 

C 

C THIS  SUBROUT  INI  CALCUIATES  THE  NOD  1 1 I!  D CHOI  I SHY  FACTORS  Of  THE 
C IIERS  URD1  R COKRM  AIIMN  NA  IRIX  MIR  A MUUINO  AVI  RAM  PROCI  SS  OF 

r r»Ri>i  r nw.  ii  ?:•  iiatb  im  iai  cui  am  imi  cokri  .vi  unimno  (himiimnis 

C AND/UR  AN  ORTMOUONAl  TRANS1  OkHAT ION  OF  INPUT  HD  MOVING  AVERAGE  DATA. 

C 


non 


piumvirm  y< irousi  > , hi  ta<nq>  , r<nii>  , uki  <nu>  . di  irons?) , 

1 Al  ( IKON.'."',  IK0NS2) ,E(  IKDNS1  > 

DA  r A (INI  /I  .0/ 

DA  I A LT'S/I.E-IO/ 

c 

C EPS  JS  CR  111  R I ON  FOR  TFSTJNG  FOR  7fl!0  DIAGONAL 

c 

C TEST  IOR  INVALID  I'AKAME  TERS 
C 

IF  AID  T-1 

II  (NU.LT.1 ) GO  TO  29? 

HAUL  r->2 

IK(  < IOPTB.NF.O)  . AND.  ( I0PTB.NE.1  > > GO  TO  299 
IFAUL  T-3 

IF(  < IOI'TE.NE.O) -AND.  dOPTE.NF.  I >)  GO  TO  299 
IFAULT-4 

IFdRONS2.LE.N0>  GO  TO  29? 

AUTOCORRELATIONS  AND  INITIALIZATION 

IF ( I OP TP . EQ . O)  CALL  NACV1 ( NO , BETA , SIG5Q , R . RO , IF  1 ) 

D< I ) - ONE 
NQP1 -NQ+1 
DO  ■S  I-1.N0 
Rd)«Rd)/RO 
3 CONTINUE 

DO  10  1-1 . NOPI 
AL  1 1 . I )— ONE 
10  CONTINUE 

IFdOPTE.E0.il  Eli  >— V ( 1 ) 

C ROW  K OF  FACTORIZATION 
C HI  IS  NUMBER  OF  LEADING  ZEROS 

C M3  IS  NUMBER  OF  ELEMENTS  TO  CALCULATE 

C MR  IS  SHIFTING  FACTOR  IN  STORAGE  OF  PREVIOUS  RONS 
C INDD+J-1  IS  RON  NUMBER  OF  AL  FOR  CALCULATING  UKL(J) 

C 


20 

30 

40 

30 


DO  130  K-2, ITERS 
M1-MAXOIK-N0-1 ,•> 

M3-MIN0IK-1 . NQ) 

INDR-K-H1-1 

MR-MAX0<K-N0-2,0> 

INDD-H1 +1 -MR 

UKL ( 1 > -R  < INDR  > /D< INDD > 

IFIM3.LT. 2)  GO  TO  4* 

DO  30  J-2.M3 
JM1-J-1 

INDR1 -INDR-JH1 
C-RdNDRI  ) 

INDLR-INDD+ JM1 
M4— MAXO(  Ml  ♦ J-NQ , 1 > 4M1  +1 
DO  20  1-1,  JM1 
INDD1 -INDD+I-1 
INDLC-M4-I 

C-C-DI INDD I )»WKL< I )*AL < INDLR , INDLC) 
CONTINUE 

UKL ( J)-C/D( INDLR) 

CONTINUE 

CONTINUE 

DK-ONE 

DO  30  J-1  ,M3 
INDD- Ml ♦J-NR 

DK -DK- P ( I NDD ) -NKL I J) *NKL I J) 

CONTINUE 

IF  1 1UPTE.EQ.G)  00  TO  AS 
C-Y(K) 

DO  AO  1-1, M3 
I NPUKL -M3- J ♦! 


' 


f r 

:■ 


hjJt 

•v. 

> 


I 

I' 


inpi:  -k-i 

C-l  WKL  d NPWNL  > 4E d NDC  > 
AO  CONTINUE 
E<K>-r: 


o n n 


A3  CONTINUE 
lFAIILT-3 

IF (DK.LT .EPS)  GO  TO  19? 

IFIK.UT.NQFI > GO  TO  BO 
D<  K > -DK 

DO  70  1-1  ,i  . 

ALIK,  I J=UKL<I> 

70  CONTINUE 

GO  TO  130 
BO  CONTINUE 

DO  82  I-l .NO 

IP1-I41 

D( I )=D< IP1 ) 

DO  81  J=1  .NO 
AL(I, J>=AL(IP1 , J> 

81  CONTINUE 

82  CONTINUE 
DINQP1 >=DK 

DO  83  1=1 .NO 
AL ( NQP1 , I)-UKL(I) 

83  CONTINUE 

DEE  IS  CONVERGENCE  CRITERION 

IFPI0PT8.E0.0)  GO  TO  110 
DO  90  1*1  , NQ 

I F < ABS  < AL  < NOP 1.1) -AL  < NQ , I ) > . GE . DEL ) GO  TO  130 
90  CONTINUE 

IF ( RO-ABS ( D ( NQP1 ) -D<  NQ ) ) . GE . DEL ) GO  TO  130 
DO  100  1=1 , NQ 
INDUKL-NQ-I-M 
BETA! I I-UKL  < INDUKL  > 

100  CONTINUE 

SIGSQ=R0»D<NQP1  > 

IFAULT-0 
GO  TO  130 


110  CONTINUE 

DO  120  1*1 , NQ 
INDB-NQ-I41 

IF<ABS<UKL<  I )-BETA(INDB)).GE.DEL)  GO  TO  130 
120  CONTINUE 

IF<ABS<D<N0P1 >»RO-SIGSQ> .LT. DEL)  00  TO  130 
130  CONTINUE 
IFAULT-A 

IF < IOPTB.EQ. 0)  GO  TO  199 
DO  140  I-1.NQ 
INDUKL -NQ-I41 
BETA< 1 >-UKL (INDUKL  > 

140  CONTINUE 

SIGSO-R0ODK 

00  TO  199 
130  CONTINUE 

IFAULT-0 

lF(IOPTE.EO.O)  GO  TO  199 
1FIK.E0. ITERS)  GO  TO  199 
KP1-K41 

DO  170  J-KP1, ITERS 
C*V< J) 

DO  1 AO  I-1.N0 
INDf-J-I 

C=C-BCTA<1>4C<INDE> 

I AO  CONTINUE 

E<  J1  =C 

170  CONTINUE 
II  AUL1-0 
1*9  CONTINUE 

DO  200  I-1.N0 
RID- RdloRO 
200  CONTINUE 
299  RtlllKN 

1 ND 

SUBROUT  I NE  H AC  V 1 < NO  . BE  T A , S I It.VQ  , R , NO , I F AUL  T ) 


o n n n 


! 

I 


f 


THIS  SUBROUTINE  CALCULATES  AUTOCOVARIANCES  OF  A HOVING  AVERAGE 
PROCESS  OF  ORDER  NQ  GIVEN  ITS  COEFFICIENTS  AND  RESIDUAL  VARIANCE. 

DIMENSION  BET  A < NQ  > , R < NQ  > 

DATA  ONE  /1.Q/ 

IFAULT-1 

IF(NQ.LT.I)  GO  TO  40 
C-ONE 

DO  10  1*1 , NQ 
C»C*BETA(I)*BETA<I> 

10  CONTINUE 
t RO*SIGSQ»C 

DO  20  IV-1 , NQ 
NQNIV-NQ-IV 
C*BETA< IV) 

DO  30  1*1 , NQHIV 
IVPI-IV*I 

C-C+BETAt I ) »BETA( IVPI > 

30  CONTINUE 

R<IV>*C«S10SQ 
20  CONTINUE 

C 

C 

IFAULT-0 
40  RETURN 
END 


'S 
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